home *** CD-ROM | disk | FTP | other *** search
/ EnigmA Amiga Run 1995 November / EnigmA AMIGA RUN 02 (1995)(G.R. Edizioni)(IT)[!][issue 1995-11][Skylink CD].iso / earcd / program / gcc / libmat.lha / src / meig.cc < prev    next >
C/C++ Source or Header  |  1980-01-01  |  327b  |  22 lines
  1. //                MATRIX LIB
  2. //            TOMMY JOHANSSON 1995
  3.  
  4. #include "matrix.h"
  5. Matrix eigen(const Matrix & A)
  6. {
  7.     Matrix C(A);
  8.     Matrix L(A.m);
  9.     Matrix R(A.m);
  10.     Matrix temp(A.m,1);
  11.     int i,ch=A.m;
  12.     while(ch>1)
  13.     {
  14.         if(fabs(C.koff[ch][1])<EPS) ch--;
  15.         LR(C,L,R);
  16.         C=R*L;
  17.     }
  18.     for(i=1;i<=A.m;i++)
  19.         temp.koff[i][1]=C.koff[i][i];
  20.     return(temp);
  21. }
  22.